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Abstract: In this paper, the enhanced (GIG )-expansion 
method is used to assemble the traveling wave solutions 
involving parameters of the Foam Drainage equation, where 

G = G(g ) satisfies a second order linear differential 

equation. The traveling wave solutions are expressed by the 
hyperbolic functions and the trigonometric functions. When 
some arbitrary functions included in these solutions are taken 
as some special functions, these solutions possess profuse 
structures. This method is direct, concise elementary and 
effective and can be used for many other nonlinear evolutions 
equations. 
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1. Introduction: 

The innovation of the soliton, its remarkable properties and the 
incredible richness of structure are all included in its 
mathematical description. The story begins with the observation 
by Jhon Scott Russell of "the great wave of translation" He 
states: 

"Its height gradually diminished and after a chase of one or two 
miles I lost it in the windings of the channel. Such, in the month 
of August 1834, was my first chance interview with that singular 
and beautiful phenomenon which I have called the Wave of 
Translation, a name which it now very generally bears [1]". 

The appearance of solitary wave solutions in nature is quite 
common. References can be made to Bell-Shaped sech- solutions 
and kink-shaped tanh-solutions model wave phenomena in 
fluids, plasmas, elastic media, electrical circuits, optical fibers, 
chemical reactions etc. The traveling wave solutions of the 
Korteweg-de Varies (KdV) and Boussinesq equations, which 
describe water waves, are famous examples as well. 

In recent years, the exact solutions of nonlinear PDEs have been 
investigated by many authors (see for examples [2-26]) who are 
interested in nonlinear physical phenomena. Many powerful 
methods have been presented such as the extended tanh-function 
method [2, 3], the tanh-function method [4, 5, 6, 7], the F- 



function expansion method [8, 9], the exp-function expansion 
method [10, 11, 12, 13], the Jacobi elliptic function expansion 
method [14, 15], the Backlund transform [16], The homogeneous 
balance method [17], the homotopy analysis method [18, 19], 

the(G7G) -expansion method [20, 21, 22, 23], the exp(-0(^))- 
expansion method [24, 25, 26] and so on. 

In this paper, we contemplate the foam drainage equation and 

apply to the enhanced ( G I G )-expansion method. Foams are of 
huge significant in many technological processes and 
applications and their properties are subjects in rigorous studies 
from both sensible and scientific point of view. Foams are 
general in foods and personal care products such as lotions and 
creams and foams over and over again arise for the period of 
clearing of clothes and scrubbing [27]. There are now a lot of 
applications of polymeric foam [28] and more newly metallic 
foams which are foams throughout metals such as aluminum 
[29]. We give attention to a quantitative explanation of the 
coupling of drainage. Foam drainage is the flow of liquid 
through plateau borders and intersections of four channels 
between the bubbles driven by gravity and capillarity. 

The paper is organized as follows. In section 2, we describe 
briefly the enhanced ( G'/G)-expansion method, where 
G = G(<^) satisfies a second order linear differential 

equation G" + AG — 0 , where ^ = k ( X + CD t ) , where 

k , CD and X are constants. The degree of this polynomial can be 

determined by considering the homogeneous balance between 
the highest order derivatives and the nonlinear terms appearing 
in the given nonlinear equations. In section 3, we apply this 
method to the foam drainage equation. In section 4, we 
discussed the physical elucidation and graphical demonstration. 
Finally, in section 5, some conclusions are given. 

2. Description of the enhanced ( G I G)-expansion method: 

In this section, we portray in niceties the enhanced (G'/G)- 
expansion method for outcome traveling wave equations of 
nonlinear equations. Any nonlinear equation in two independent 
variables x and t can be articulated in the following form: 
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M(u,u t ,u x ,u tt ,u xx ,u xt ) =0 



(2.1) 



where u(Jg) = u(x,t) is an anonymous function, 9? is a 

polynomial of u(x,t) and its partial derivatives wherein the 

uppermost order derivatives and non-linear stipulations are 

concerned. The following steps are concerned in result the 
solution of nonlinear Equation (2.1) using this method. 

Step 1: The given PDE (2.1) can be distorted into ODE using the 
transformation lg — X + COt , where CO is the speed of traveling 

wave such that CO E R — {0} . 

The traveling wave alteration permits us to reduce Eq. (2.1) to 
the following ODE: 



*¥{u,U ,U , )= 0, 



(2.2) 



Where ^is a polynomial in u(jg) and its derivatives, where 

, du „ d u 
u (g) , u (g) - ——^ , and so on. 



Step 2: Now, we expect that the Eq. (2.2) has a general solution 
of the form 



(2.3) 



(0 = 2 



i= n 



a^G'/Gj 



-^{g'/gJ- 1 . 



(i+a{g'/g)) 

G = G(^) satisfy the equation 



i+ 



(G7G) 



2\ 



subject matter to 
the condition that 



(2.4) 



Where a i ,b i {-n <i <n\n e N), and X are constant to be 
determined, provided that c = ±land jU ^ 0 . 

Step 3: The positive integer n can be determined by balancing 
the highest order derivatives to the highest order nonlinear terms 
appear in Eq. (2.1) or in Eq. (2.2). More accurately, we define 

the degree of u{jg) as D(u(<^)) = n which gives respond to 

the degree of other expression as follows: 



D 



( d q u^ 


= n + q, D 


f 


d q u 






u p 












J 



= np + s{n + q), (2.5) 



Step 4: We substitute Eq. (2.3) into Eq. (2.2) and use Eq. (2.4). 
We then collect all the coefficient of (G'/G) j and 



(G7G) 



V 



! , (G'lG) 



2 \ 



together. Since Eq. (2.3) is a 



J 



solution of Eq. (2.2). we can set each of the coefficient equal to 
zero which leads to a system of algebraic equations in terms of 
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a i ,b i {—n < i < n\n e N)> A and CO . One can solves easily 
these system equations using Maple. 

Step 5: For jU < 0 general solution of Eq. (2.4) gives 



and 



G 
G' 



= J^ju tanh(A + -J^ju 
= V-// coth(A + 7-// 



(2.6) 



(2.7) 



And for // > 0 , we get 



<7 
G 



= Jjutan(A-y[ji £), 



(2.8) 



and 



G 



= y[JlCOl(A + JJl £), 



(2.9) 



where A is an arbitrary constant. Finally we can construct a 
number of families of travelling wave solutions of Eq. (2.1) by 

substituting the values of a t , b t {—n <i<n;n<= N), X and 

CO (obtained in Step 3) and using Eq. (2.6) to Eq. (2.9) into Eq. 
(2.3). 

3. Application: 

In this section, we pertain the enhanced ( G' I G )-expansion 
method to put up the traveling wave solutions for foam drainage 
equation 



(3.1) 




Where y/ is the cross segment of a channel formed where three 

films meet, usually indicate as plateau border and X and t are 
scaled position and time coordinates respectively. By means of 

the traveling wave variable I// = (p(g ) , g = k (x + COt) , where 

k and CO are constant, Eq. (3.1) is altered into the following 
ordinary differential equation (ODE) for U — U (g) . 



k r- 

kcocp' + k(cp 2 — y<p<p'y = 0, 

2> 



(3.2) 



Where prime denote the differentiation with regard to g . Now, 

integrating Eq. (3.2) with respect to<^ and bearing in mind the 
zero constant for integrating we have, 



k i— 

kcocp + k{(p 2 — ^cpcp') = 0, 

2> 



(3.3) 
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Afterwards we use the transformation 

(p{$) = u 2 ($) , 

to exchange Eq.(3.3) to 

kcou 1 +ku A — k 2 u 2 u = 0 , 

or equivalently 

co + u 2 —ku — 0, 



(3.4) 



(3.5) 



(3.6) 



Now, balancing the highest order derivative of U ' and nonlinear 

2 1 

term of U in Eq.(3.6), we obtain Ti — l. 

So the equation (3.6) takes the following solution 



u (£) = a Q + 



\+x{g'ig) (g'/g) ov 7 | [ // 



+2>i Jo - 



u {g'/g) 



2\ I / 

+ ^(G7G)- 2 ja 



1+ (g/g) 



2\ 



(3.7) 



where G = G(<^) satisfy the Eq. (2.4). Substituting Eq. (3.7) 
hooked on the Eq. (3.6) and using Eq. (2.4), we get a polynomial 



in (G7G) J and(G7G) 7 \CJ 



V 



! , (G'/G) 



2 \ 



J 



Surroundings the coefficient of 



(g'/g) 



and 



(G7G> 



a 



v 



1+ (G7G) 



2 A 



equal to zero, we obtain a 



system containing a bulky number of algebraic equations in 
terms of unknown coefficients. We have solved this system of 
equations using Maple 13 and obtained the following set of 
solutions: 



Set 



co-k ju,a 0 =-k fiA, a x = 0, a x = kju,b Q =0,b l = 0, b A - 0 



Set 



2 2 

a = k [i,a 0 =k fiX, flj = -k fii -k, a A = 0, b 0 = 0, b x = 0, b A = 0 



Set 



1 



1 



1 



1 ftk 



a=-k //, ^^--kfil, a { = 0, a_ { =-kfi, \ =+-—, \ = 0,b_ l =0, 



4 



2 



2 



2 V<7 



Substituting Set 1-Set 3 into Eq. (3.7) along with Eq. (2.6)-Eq. 
(2.9) and using the transformation \j/{x, t) = (p{%) = u (g) , 
we have the following families of traveling wave solutions: 
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Hyperbolic function solutions: 



when jU < 0 , we get the following three families of hyperbolic 
function solutions. 

J*— c oth( A + yf^ju %) 



Family 1 : y/ l (x, t) = 



V 



J 



V 



J^JL tanh (A + yf^ju E) 



where ^ = k(x-\-k /Jt), 

f I / ^2 

k juA — k yj— ju tanh( A + v ~ M <?) 



Family 2: y/ 3 (x,t) = 



V 



J 



f 



y/ 4 (x,t) = 



1 + /l^/- ju tanh(A + J- ju £) 



k juA-k ^- ju coth( A + ^- ju 



V 



1 + A 7-// coth( A + 



y 



where E, - k (x + k /Jt), 



Family 3: y/ 5 6 (x,t) = 



k ^ (coth( A + 7^ f) ± esc h(A + £)) 



V 



2/ 



J 



2/ 



(tanh( A + 7^// ^ ) ± / sec h(A + J^/z ^)) 



\ 2 



1-2 



where <^ = /: (x + — /: /i f ) 



Trigonometric function solutions: 



when /J > 0 , we get the following three families of 
trigonometric function solutions. 

Family 4: l// 9 (x, t) = {k y[JlCOt(A - yfju 

y/ m (x, t) = [ky[Jt tan(A -Jjug) 



1 2 

i - : where <^ = k (x + /Jt), 



Family 5: l// n (x,t) = 



V 



k juX-k^ 


tan (A - 






1 + A/ 


tan (A-yj 





J 



¥\2 (*> t) = 



( 



V 



k juA-kyj 


^//cot(A + ^/ 




1 + /L^/ 


^//cot(A + ^/ 





where <^ = k (x + jut), 



f 



Family 6: y/ uu (x,t) = 



kyJM 



(cot (A-y[ju^)± csc (A + Jju € )) 



V 



^15,16(^0 = 



V 



2 



(tan (A + 7// + sec ( A + ^ f)) 
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where E — k (x + — k at) , 

4 

Remark: All the obtained solutions have been tartan with maple 
by putting them back into the original equations and found 
correct. 

4. Results and Discussion: 

In this section, we will portray the physical elucidation and 
graphical demonstration of the solutions of the Foam Drainage 
equation. 

4.1. Physical Elucidation: 

A quantum of energy or quasi particle that can be propagate as a 
traveling wave in nonlinear systems and is neither preceded nor 
followed by another such disturbance; does not obey the 
superposition principle and does not dissipate. Soliton waves can 
travel long distances with little loss of energy or structure. 
Solitary waves are localized travelling waves travelling with 
constant speeds and shape, asymptotically zero at large 
distances. Solitons are special kinds of solitary waves. The 
soliton solution is especially localized solution, hence 

u (<!;), u"(<!;) and u"\^) — Q as with^ = ±00,^ = x — ct . 

Solitons have a remarkable soliton property in that it keeps its 
identity upon interchanging with other solutions. Solitons 
solutions are give rise to particle like structures, such as 
magnetic monopoles etc. So, soliton are everywhere in the 

nature. The solution l// 1 (x,t) is represented the soliton 

solution. The Fig. 1 has been shown the shape of the solution 

y/ x {x,t) for ju = —l,k = 0.5, A = 0 

with — 10 < x,t < 10. For the fixed values 
/J = — 1, k = 1, A — 0 with — 3 < X, t < 3 , solution 

y/ 2 (x, t) represents the exact solitary wave solutions of bell 
type which is shown in Fig. 2. For the fixed values 
ju = 1, k = 1, A = 0 with - 3 < x, t < 3 , solution y/ l0 (x, t) 

represents periodic wave solutions. Periodic solutions are 
traveling wave solutions are periodic such as COS(x — t) . The 
exact periodic traveling wave solution is represented graphically 
in Fig. 3. Finally, solution y/ 5 (x,t) is singular soliton of the 
Foam Drianage equation. The Fig. 4 shows the shape of the 
singular soliton solution, obtain from l// 5 (x,t) for 

ju = -1, k = 1, A = 0 with - 3 < x, t < 3 . 

4.2. Graphical Demonstration: 

In this sub-section, the graphical demonstration of the solutions 
is given below in the figures (Fig. 1-4) with the aid of 
mathematical software Maple 13. 



xlO 8 - 



xlO 8 - 



xlO 8 - 



xlO 8 - 




-100 



-31- 



1 1 j 1 1 1 1 1 1 

5 10 



Fig. 1. Soliton profile of Foam Drianage equation for 
JU = -1, £ = 0.5, A = 0 with- 10 < x,t <10. (Only 

shows the shape of y/ x (x, t) ), The left figure shows the 3D plot 
and the right figure shows the 2D plot for t — 0 . 
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Fig. 2. Bell type soliton profile of Foam Drainage equation for 
/J = —1, k = 1, A = 0 with — 3 < X, t < 3 . (Only shows the 

shape of y/ 2 (x, t) ), The left figure shows the 3D plot and the 
right figure shows the 2D plot for t — 0 . 
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Fig. 3: Periodic type soliton profile of Foam Drianage 
equation for jU = 1, k = 1, A — 0 with — 3 < X, t < 3 . 

(Only shows the shape of lf/ l0 (x, t) ), The left figure shows 

the 3D plot and the right figure shows the 2D plot for t = 0 . 
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Fig. 4: Singular soliton profile of Foam Drianage equation for 
/J = —1, k = 1, A = 0 with — 3 < X, t < 3 . (Only shows 

the shape of y/ 5 (x, t) ), The left figure shows the 3D plot and 
the right figure shows the 2D plot for t — 0 . 

5. Conclusion: 

Foaming occurs in many cleansing and amalgamation processes. 
The drainage of fluid foams involves the dealings of gravity, 
surface tension and viscous forces. In this paper, we use the 

enhanced ( G I G )-expansion method to seek the traveling wave 

solution of the foam drainage equation. The enhanced ( G I G )- 
expansion method provides a very effective and powerful 
mathematical tool for solving nonlinear equations in 
mathematical physics. 
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